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Abstract. Muon spin spectroscopy and in particular the avoided level crossing technique 
is introduced, with the aim of showing it as a very sensitive local probe for electron spin 
relaxation in organic semiconductors. Avoided level crossing data on TMS-pentacene at different 
temperatures are presented, and they are analysed to extract the electron spin relaxation rate, 
that is shown to increase on increasing the temperature from 0.02 MHz to 0.33 MHz at 3 K and 
300 K respectively. 



1. Introduction 

Organic semiconductors are extremely promising materials for spintronics thanks to their long 
spin coherence times. However, spin phenomena in these materials are not yet fully understood, 
partly due to the lack of suitable characterisation techniques, since those tuned for inorganic 
semiconductors are not always successful for organic materials [HE]- Muon spin relaxation (/uSR) 
has been proven to be an applicable probe in this field, as it is sensitive to spin dynamics in 
organic semiconductors on a molecular lengthscale [3] . In organic semiconductors the implanted 
muons form a hydrogen-like system, called muonium, by capturing an electron. In muonium the 
muon and electron spins are coupled through hyperfine interaction. As a consequence relaxation 
on the electron spin can induce changes in the muon spin, which can be detected. Here we show 
how muon spin relaxation and in particular the avoided level crossing (ALC) technique can be 
used to probe the electron spin dynamics in organic semiconductors. 

2. Muon spin relaxation and avoided level crossings 

Muon spin relaxation is an experimental technique using positive muons as microscopic spin 
probe. Positive muon has a finite lifetime r ~2.2 ^s, and it radioactively decays emitting a 
positron and two neutrinos according to the process /x"*" — )• e+ + z/g + z^. The positron is emitted 
preferentially in the direction of the muon's spin at the time of the decay, due to parity violation 
in weak interaction [H [5] . As a consequence measuring the time dependence of the emission 



direction of the positron provides information on the time evolution of the muon spin. 

In ;uSR experiments 100% spin polarised muons are implanted in the sample. Once implanted 

in semiconductors the muon gives rise to a series of ionisation and scattering processes, through 

which its energy decreases from 4MeV to a few keV. Then it takes part in a series of electron 

capture and loss processes that reduce its energy further. Finally it thermalises to form 

interstitial positive muon or it captures an electron and form a hydrogen-like species called 

muonium (Mu) [5j. Thermalisation of muons is fast (<ns) compared to its lifetime. 

If an external magnetic field is applied in the direction of the initial muon spin, as it is in so 

called longitudinal field experiments, the Hamiltonian for muonium can be written as 

H = --i^ni-B + -fefiS-B + hS-A-I (1) 

where 7^ and 7e are the muon and electron gyromagnetic ratios, B is the applied longitudinal 
field, I and S are the muons and electrons spins and A is the hyperfine tensor [5]. The energy 
levels as a function of the applied magnetic field that originates from this Hamiltonian for an 
isotropic hyperfine interaction are shown in the Breit-Rabi digram in Figure ma). At B=0 a 
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Figure 1. (a) The muonium Breit-Rabi diagram with a hyperfine coupling constant of A=80 
MHz (b) Zoom-in region of the energy levels near the ALC. 



singlet and a triplet are obtained, the degeneracy of the latter being lifted by the magnetic field. 
At high enough fields (Paschen-Back regime, where most of the experiments are performed) the 
electron and muon spins are decoupled and the state of the system can be described by products 
of pure Zeeman states for the electron and muon. 

At a high field which depends on the value of the hyperfine coupling constant, two of the levels 
should cross, as can be observed in Figure[T]^b) (dashed lines). However, if further interactions are 
present, as an example with the nuclei, this crossing can be avoided. In fact, for a Hamiltonian 
H = Hq + H^ where can be treated as a perturbation, the correction to the energies up to 
second order is given by 

m n 

where and Xn are the n-th zero order energy and eigenstate and the denominator becomes 
tiny when the two unperturbed levels approach The effect of this correction to the energy 
levels is to avoid the crossing (See Figure [T]^b)). Close to this avoided level crossing the states are 
no longer pure Zeeman states, but they are instead mixtures of states corresponding to different 
values of the muon magnetic quantum number. As a consequence a loss of polarisation of the 
muon can be observed [U [5] . ALCs can be classified by the change in total quantum number 
M. If the hyperfine interaction is isotropic, the most common kind of ALC is the one with 



AM=0, that corresponds to a muon-nuclear spin flip-flop. However in solids the most intense 
crossing is AM=1, corresponding to a muon spin flip and is observed only when anisotropy in 
the interaction is present. AM=2 transitions exist, but they are usually very weak and narrow 
[3] . A simulated ALC curve is shown in Figure [2] (black line) . 




Figure 2. Simulated ALC resonance for different electron spin relaxation rates (eSR) 



As can be observed, the polarisation is plotted versus the applied magnetic field. The polarisation 
exhibits a deepening, representing the loss of polarisation around the ALC field. The shape and 
position of this ALC curve depends on the hyperfine coupling constants (both isotropic and 
anisotropic). However, the presence of an electron spin relaxation rate (eSR) also has an effect 
on the ALC curve. ALCs simulated using the same parameters with an increasing electron 
spin relaxation rate are shown in Figure [2j It is clear that at low electron spin relaxation rates 
neither the position nor the shape of these ALCs show substantial changes, but the amplitude 
enormously increases. 

This can be modeled through a density matrix formalism, that is also used to simulate the 
experimental data and extract the eSR. A density matrix is defined as 

p{t) = \[l+ Pit) . a + pe{t) -T + Y. P''{t)a^r'^] (3) 

where p{t), Pe{t) and P^^{t) are the muon, electron and mixed polarisations, respectively. The 
time evolution of p{t) is described by the equation 

.Sf = (4) 

from which the time dependent polarisations can be obtained [6J . The general solution can only 
be determined numerically. If an electron spin relaxation mechanism is present in the system 
its effect can be taken into account phenomenologically by introducing extra terms —\pe and 
—Xpj^ in the differential equations for the polarisations obtained from equation |4] [7]: 

t^ = ...-Ape,^ = ...-AP^-^ (5) 

dt ^ dt ^ ' 



3. Experimental methods 

(Trimethyl)silyl-pentacene (TMS-pentacene) was synthesised according to a published procedure 
[SI [2] and was further purified by repeated recrystalhsation from dicliloromethano/ethanol. Its 
structure is shown in Figure [3j The presence of the side groups on the central spine prevents the 




Figure 3. Molecular structure of TMS-pentacene. Hydrogen has been omitted for clarity. 



molecule to take the usual herringbone packing of pentacene; as a consequence the molecules 
have a face-to-face packing that allows better vr-orbitals overlap and consequently better charge 
carrier mobility [8l[9]. The measured sample was polycrystalline powder. 

Muon experiments were performed at the HIFI spectrometer at ISIS, Rutherford Appleton 
Laboratory. About 250 mg of TMS-pentacene was placed in an envelope (17mm x 17mm) 
made of 25 /um thick silver foil (99.99 % pure). Another layer of silver foil was placed in front 
of the sample as a degrader, to assure that the muons stop in the sample. 

The software Wimda [lOj was used to extract the time-integrated data and plot them as a 
function of applied field. ALCs were obtained in this way, and they were then modelled using 
the software Quantum based on the density matrix formalism described in the previous section 
|llj . The ALC simulations presented in the following are performed considering a muon-electron 
system. No nuclei are taken into account as their infiuence was discovered to be negligible. 
Monte-Carlo simulations with powder averaging and 10000 iterations for each magnetic field 
was used. 

4. Results and discussion 

Figure |4] shows the ALCs measured for TMS-pentacene at 3K, 150K and 300K. Data were also 
acquired at lOK and 25K (data not shown). ALCs are observed around 0.28T, corresponding 
to a hyperfine coupling of around 70MHz. The position and width of the ALCs do not show a 
remarcable temperature depence, while the amplitude considerably increases on increasing the 
temperature. From the time-differential data and the field dependent relaxation rate (data not 
shown, see ref [12]) this increased ALC depth is clearly related to a higher relaxation rate of the 
muon spin. 

The measured ALCs look as the superposition of three ALCs very close to each other, the 
difference arising from slightly different possible bonding angles of the muon in the molecule. As 
a consequence three contributions were also considered in the modeling. Their relative amplitude 
is determined from the 3K data, setting the eSR to a very low value. This assumption is justified 
if the electron spin relaxation is driven by an effect strongly dependent on temperature, as a 
coupling to vibrations [12]. The result is shown in Figure |4] (black line). The obtained relative 
amplitudes are 0.17, 0.24 and 0.59. These parameters are used to scale the three contributions 
at higher temperatures. An eSR different from zero is necessary in the model to account for 
the higher intensity of the ALC at higher temperatures, whereas there is only a small change 
in the hyperfine coupling constants. The results of these simulations are shown as solid lines in 




Figure 4. Muon data around the ALC. Black lines are the result of modelling. 



Table 1. Simulation parameters for TMS-pentacene. 





3K 






lOK 








ALCl 


ALC2 


ALC3 


ALCl 


ALC2 


ALC3 


A (MHz) 


61.7±0.2 


71.2±0.2 


81.9±0.2 


61.7±0.2 


71.2±0.2 


81.9±0.2 


Dl (MHz) 


3.6±0.3 


6.0±0.4 


5.8±0.4 


3.6±0.3 


6±0.4 


5.0±0.4 


D2 (MHz) 


3.6±0.5 


3.0±0.4 


4.0±0.4 


3.6±0.5 


3.0±0.5 


4.0±0.5 


eSR (MHz) 


0.02±0.02 


0.02±0.02 


0.02±0.02 


0.04±0.02 


0.04±0.02 


0.04±0.02 




25K 






150K 








ALCl 


ALC2 


ALC3 


ALCl 


ALC2 


ALC3 


A (MHz) 


62.5±0.2 


71.8±0.2 


81.6±0.2 


63.0±0.2 


72.2±0.2 


81.0±0.2 


Dl (MHz) 


3.6±0.3 


4.5±0.4 


5.0±0.4 


3.6±0.3 


4.5±0.5 


6.0±0.4 


D2 (MHz) 


3.6±0.5 


2.0±0.5 


4.0±0.5 


3.6±0.5 


2.0±0.5 


4.0±0.5 


eSR (MHz) 


0.11±0.02 


0.11±0.02 


0.11±0.02 


0.13±0.02 


0.13±0.02 


0.15±0.02 


300K 




ALCl 


ALC2 


ALC3 








A(MHz) 


66.0±0.2 


73.0±0.2 


79.0±0.2 








Dl(MHz) 


3.6±0.3 


4.5±0.4 


6.0±0.4 








D2(MHz) 


3.6±0.5 


2.0±0.5 


4.0±0.5 








eSR(MHz) 


0.33±0.02 


0.33±0.02 


0.33±0.02 









Figure |4j Values for hyperfine coupling constants A, Dl and D2 ^llj and for eSR as obtained 
from these simulations are reported in Table [T] The three ALCs tend to merge on increasing the 
temperature. These small changes in the values of A, Dl and D2, of the order of few percent, are 
likely to be due to changes in bond length at different temperatures. On the other hand there is 
a substantial increase in the eSR as a function of temperature, as can be seen in Table [T] and in 
Figure [5| Figure [5] shows the electron spin relaxation obtained from the simulation as a function 
of inverse temperature. This behaviour closely resembles the one observed in other organic 
semiconductors, like rubrene, Gaqs and TIPS-pentacene p!2|. The electron spin relaxation rate 
in rubrene and TIPS-pentacene exhibits an Arrhenius-like temperature dependence, with two 



characteristic energy scales comparable to the energy of the vibrations in these systems [1 2j . 
Data presented here are not detailed enough to extract a full temperature dependence of the 
eSR, but they are very consistent with [12j. It is also clear that the ALC spectroscopy is a 
very sensitive technique for the study of the electron spin relaxation in organic semiconductors. 
Further experiments are needed to get a detailed temperature dependence of eSR in TMS- 
pentacene. 
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Figure 5. Temperature dependence of the electron spin relaxation rate obtained from the 
simulations of the ALCs. 



5. Conclusions 

The present paper reviews the basic physics behind the avoided level crossing technique. It 
was shown that muons are a very effective probe to measure electron spin relaxation in organic 
semiconductors thanks to their coupling to the electron via hyperfine interaction. The electron 
spin relaxation rate can in fact be extracted from ALC measurements. As an example data on 
TMS-pentacene at three different temperatures are shown and analysed to obtain the electron 
spin relaxation rate, that is shown to increase on increasing the temperature. 
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